home *** CD-ROM | disk | FTP | other *** search
/ Aminet 40 / Aminet 40 (2000)(Schatztruhe)[!][Dec 2000].iso / Aminet / mus / misc / DelfScope.lha / DelfScope / src / scope_dsp.a56 < prev    next >
Encoding:
Text File  |  2000-10-03  |  33.3 KB  |  1,035 lines
  1. ;*****************************************************************************
  2. ;
  3. ;    DelfScope - oscilloscope/analyzer for Delfina DSP
  4. ;    Copyright (C) 2000  Michael Henke
  5. ;
  6. ;    This program is free software; you can redistribute it and/or modify
  7. ;    it under the terms of the GNU General Public License as published by
  8. ;    the Free Software Foundation; either version 2 of the License, or
  9. ;    (at your option) any later version.
  10. ;
  11. ;    This program is distributed in the hope that it will be useful,
  12. ;    but WITHOUT ANY WARRANTY; without even the implied warranty of
  13. ;    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  14. ;    GNU General Public License for more details.
  15. ;
  16. ;    You should have received a copy of the GNU General Public License
  17. ;    along with this program; if not, write to the Free Software
  18. ;    Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  19. ;
  20. ;*****************************************************************************
  21.  
  22.  
  23.  
  24.  
  25.     include 'delfequ.inc'
  26.  
  27.     org    p:
  28.  
  29. ;->prog+0
  30.     jmp    convert
  31.  
  32.  
  33.  
  34. ;->prog+2
  35.  
  36. ;***************************************************************************
  37. ;*** DelfModule - receives the input data
  38.     move    y:(r0),r0        ;input pipe
  39.     move    y:(r1),r1        ;output pipe
  40.     move    y:y_bufpt,r4
  41.     movec    #2048-1,m4        ;m4 - modulo addressing
  42.     do    #128,mod_receive_loop
  43.     move    l:(r0)+,x
  44.     move    x,l:(r4)+
  45.     move    x,l:(r1)+
  46. mod_receive_loop
  47.     move    r4,y:y_bufpt
  48.  
  49.     move    y:y_mod_count,a
  50.     move    #>1,x0
  51.     add    x0,a    a,x1
  52.     teq    x1,a            ;do not count while (mod_count==-1)
  53.     move    a,y:y_mod_count
  54.     rts
  55.  
  56.  
  57. ;***************************************************************************
  58. ;*** convert PCM samples
  59. ;
  60. ;  -calculate average amplitudes (into x:avg_l, x:avg_r)
  61. ;  -plot_type==1: 24bit->8bit, half sampling rate (into l:buf8)
  62. ;  -plot_type==2: spectrum analyzer
  63. ;
  64. ; x1 ... number_of_samples
  65. ; x0 ... plot_type  (1 or 2)
  66. ; y1 ... /
  67. ; y0 ... /
  68. convert
  69.     move    x1,n0            ;numsamp            !! n0
  70.     move    x0,y:y_plot_type
  71.     move    y:y_bufpt,r0
  72.     movec    #2048-1,m0        ;m0 - modulo-addressing
  73.     move    r0,x:x_busy        ;set BUSY
  74.     move    (r0)-n0
  75.  
  76. ;******    calculate average amplitude
  77.     move    r0,r6
  78.     movec    m0,m6            ;m6 - modulo-addressing
  79.     move    #>$002000,y1        ;shift right by 10 bits
  80.     clr    b    #<0,x0
  81.     move        x:(r6),a    b,y0
  82.     do    n0,avg_amp_loop
  83.     abs    a    y:(r6)+,b
  84.     abs    b    a,x1
  85.     mpy    x1,y1,a    b,x1
  86.     mpy    x1,y1,b
  87.     add    x0,a
  88.     add    y0,b    a,x0
  89.     move        x:(r6),a    b,y0
  90. avg_amp_loop
  91.     move    y:y_plot_type,n1    ;jmp
  92.     move    #>y_plot_jmp,r1        ;jmp
  93.     move    x0,x:x_avg_l        ;must be divided by numsamp for avg
  94.     move    y0,x:x_avg_r        ;must be divided by numsamp for avg
  95.     move    y:(r1+n1),r1        ;jmp
  96.     movec    m7,m6            ;m6 - linear addressing
  97.     jmp    (r1)
  98.  
  99.  
  100. ;******    plot one: convert 24bit samples to 8bit (into l:buf8)
  101. plot_1    move    n0,a            ;numsamp
  102.     move    #>5,x0
  103.     add    x0,a    #>6,x0
  104.     move    a,x1
  105.     jsr    <delf_div        ;6 samples per loop
  106.     rnd    a    #>l_buf8,r1
  107.     ori    #$04,mr            ;SCALE DOWN mode
  108.     do    a1,p248_loop
  109.     move        l:(r0)+,ab
  110.     move        x:(r0),x1
  111.     add    x1,a    y:(r0)+,x1
  112.     add    x1,b    #>$ff0000,y1
  113.     asr    a
  114.     asr    b
  115.     and    y1,a
  116.     and    y1,b    a1,x0
  117.     move        b1,y0
  118.     move        l:(r0)+,ab
  119.     move        x:(r0),x1
  120.     add    x1,a    y:(r0)+,x1
  121.     add    x1,b    a,x1    ;scale down
  122.     mpy    x1,#8,a
  123.     move        #>$00ff00,y1
  124.     and    y1,a    b,x1    ;scale down
  125.     mpy    x1,#8,b
  126.     or    x0,a
  127.     and    y1,b    a1,x0
  128.     or    y0,b
  129.     move        b1,y0
  130.     move        l:(r0)+,ab
  131.     move        x:(r0),x1
  132.     add    x1,a    y:(r0)+,x1
  133.     add    x1,b    a,x1    ;scale down
  134.     mpy    x1,#16,a
  135.     move        #>$0000ff,y1
  136.     and    y1,a    b,x1    ;scale down
  137.     mpy    x1,#16,b
  138.     or    x0,a
  139.     and    y1,b    a1,x:(r1)
  140.     or    y0,b
  141.     move        b1,y:(r1)+
  142. p248_loop
  143.     andi    #$f3,mr            ;NO SCALING mode
  144.     movec   m7,m0            ;m0 - linear addressing
  145.     move    #<0,x0
  146.     move    x0,x:x_busy        ;clear BUSY
  147.     rts
  148.  
  149.  
  150.  
  151. ;******    plot two: spectrum analyzer
  152. plot_2    move    #>l_fft_in,r7
  153.     move    y:y_bufpt,r0
  154.     move    #1024,n0
  155.     move    #$002000,x1        ;#$002000 - shift right by 10 bits
  156.     move    (r0)-n0
  157.     move    r7,r1            ;#>l_fft_in
  158.     move    #>y_fft_window,r2
  159.     clr    b        x:(r0)+,x0
  160.     do    #1024,fft_input_left
  161.     mpyr    x0,x1,a
  162.     move    a,x0        y:(r2)+,y0
  163.     mpyr    x0,y0,a        x:(r0)+,x0
  164.     move            ab,l:(r1)+    ;set re, clear im
  165. fft_input_left
  166.  
  167. ;FFTR2CN
  168. ;
  169. ; This program originally available on the Motorola DSP bulletin board.
  170. ; It is provided under a DISCLAMER OF WARRANTY available from
  171. ; Motorola DSP Operation, 6501 Wm. Cannon Drive W., Austin, Tx., 78735.
  172. ; 
  173. ; Radix 2, In-Place, Decimation-In-Time FFT (fast).
  174. ; 
  175.     move    #>l_fft_sincos,n7
  176.  
  177.     move    r7,r0            ;#>l_fft_in
  178.     move    #1024/4,n0
  179.     move    n0,n4
  180.     move    n0,n6
  181.     movec    #1024-1,m0        ;m0 - modulo addressing
  182.     movec    m0,m1            ;m1 - modulo addressing
  183.     movec    m0,m4            ;m4 - modulo addressing
  184.     movec    m0,m5            ;m5 - modulo addressing
  185.     move            x:(r0)+n0,x0
  186.     tfr    x0,a        x:(r0)+n0,y1   
  187.     do    n0,fft_twopass
  188.     tfr    y1,b        x:(r0)+n0,y0
  189.     add    y0,a        x:(r0),x1
  190.     add    x1,b        r0,r4
  191.     add    a,b        (r0)+n0
  192.     subl    b,a        b,x:(r0)+n0
  193.     tfr    x0,a        a,x0        y:(r0),b
  194.     sub    y0,a                y:(r4)+n4,y0
  195.     sub    y0,b        x0,x:(r0)
  196.     add    a,b                y:(r0)+n0,x0
  197.     subl    b,a        b,x:(r0)
  198.     tfr    x0,a        a,x0        y:(r4),b
  199.     add    y0,a                y:(r0)+n0,y0
  200.     add    y0,b        x0,x:(r0)+n0
  201.     add    b,a                y:(r0)+,x0
  202.     subl    a,b                a,y:(r4)+n4
  203.     tfr    x0,a                b,y:(r4)+n4
  204.     sub    y0,a        x1,b
  205.     sub    y1,b        x:(r0)+n0,x0
  206.     add    a,b        x:(r0)+n0,y1
  207.     subl    b,a                b,y:(r4)+n4
  208.     tfr    x0,a                a,y:(r4)+
  209. fft_twopass
  210.  
  211.     move    #1024/8,n1        ;butterflies per group
  212.     move    #<4,n2            ;groups per pass
  213.     movec    #<0,m6            ;m6 - bit-reversed addressing
  214.     do    #10-3,fft_pass
  215.     move    r7,r0            ;#>l_fft_in
  216.     move    r0,r1
  217.     move    n1,r2
  218.     move    r0,r4
  219.     move    (r1)+n1
  220.     move    r1,r5
  221.     move    n7,r6            ;#>l_fft_sincos
  222.     lua    (r2)+,n0
  223.     move    n0,n4
  224.     move    n0,n5
  225.     move    (r2)-            ;butterfly loop count
  226.     move            x:(r1),x1    y:(r6),y0
  227.     move            x:(r6)+n6,x0    y:(r0),b
  228.     mac    x1,y0,b                y:(r1)+,y1
  229.     macr    -x0,y1,b            y:(r0),a
  230.     do    n2,fft_grp
  231.     do    r2,fft_bfy
  232.     subl    b,a        x:(r0),b    b,y:(r4)
  233.     mac    -x1,x0,b    x:(r0)+,a    a,y:(r5)
  234.     macr    -y1,y0,b    x:(r1),x1
  235.     subl    b,a        b,x:(r4)+    y:(r0),b
  236.     mac    x1,y0,b                y:(r1)+,y1
  237.     macr    -x0,y1,b    a,x:(r5)+    y:(r0),a
  238. fft_bfy    move    (r1)+n1
  239.     subl    b,a        x:(r0),b    b,y:(r4)
  240.     mac    -x1,x0,b    x:(r0)+n0,a    a,y:(r5)
  241.     macr    -y1,y0,b    x:(r1),x1    y:(r6),y0
  242.     subl    b,a        b,x:(r4)+n4    y:(r0),b
  243.     mac    x1,y0,b        x:(r6)+n6,x0    y:(r1)+,y1
  244.     macr    -x0,y1,b    a,x:(r5)+n5    y:(r0),a
  245. fft_grp    move    n1,b
  246.     lsr    b    n2,a    ;divide butterflies per group by two
  247.     lsl    a    b,n1    ;multiply groups per pass by two
  248.     move    a,n2
  249. fft_pass
  250.  
  251.     move    #<2,n0
  252.     move    n0,n1
  253.     move    #1024/4,n4
  254.     move    n4,n5
  255.     move    r7,r0            ;#>l_fft_in
  256.     move    #>l_fft_out,r4
  257.     move    r4,r2
  258.     lua    (r0)+,r1
  259.     lua    (r2)+n2,r5
  260.     movec    #<0,m4            ;m4 - bit-reversed addressing
  261.     movec    m4,m5            ;m5 - bit-reversed addressing
  262.     move    n7,r6            ;#>l_fft_sincos
  263.     move    (r5)-n5
  264.     move            x:(r1),x1    y:(r6),y0
  265.     move            x:(r5),a    y:(r0),b
  266.  
  267.     do    n2,fft_lastpass
  268.     mac    x1,y0,b        x:(r6)+n6,x0    y:(r1)+n1,y1
  269.     macr    -x0,y1,b    a,x:(r5)+n5    y:(r0),a
  270.     subl    b,a        x:(r0),b    b,y:(r4)
  271.     mac    -x1,x0,b    x:(r0)+n0,a    a,y:(r5)
  272.     macr    -y1,y0,b    x:(r1),x1    y:(r6),y0
  273.     subl    b,a        b,x:(r4)+n4    y:(r0),b
  274. fft_lastpass
  275.     move            a,x:(r5)+n5
  276.  
  277.     movec   m7,m0            ;m0 - linear addressing
  278.     movec   m7,m1            ;m1 - linear addressing
  279.     movec   m7,m4            ;m4 - linear addressing
  280.     movec   m7,m5            ;m5 - linear addressing
  281.     movec   m7,m6            ;m6 - linear addressing
  282. ;------------------------------------------------------------------------
  283.  
  284.  
  285.  
  286.  
  287. ;convert complex output values to magnitude (dB-scale)
  288.     move    #>l_fft_out,r0
  289.     move    #>l_buf8,r1
  290.  
  291.     do    #416,fft_magnitude
  292.     move            x:(r0),x0
  293.     mpy    x0,x0,a        y:(r0)+,y0
  294.     mac    y0,y0,a        #<$40,x0    ;a=yr**2+yi**2
  295.     move    a1,y1                ;transfer sum of squares to y
  296.     move    a0,y0
  297.  
  298.     clr    b        x0,x1    ;init root
  299.     do    #23,fft_mag_sqrt
  300.     mpy    -x0,x0,a        ;square and negate the guess
  301.     add    y,a            ;compare to double precision input
  302.     tge    x0,b            ;update root if input >= guess
  303.     tfr    x1,a            ;get bit to test
  304.     asr    a            ;shift to next bit to test
  305.     add    b,a        a,x1    ;form new guess
  306.     move    a,x0            ;save new guess
  307. fft_mag_sqrt
  308.  
  309.     move    #>y_fft_log_coef,r2
  310.     move    #>$000008,a
  311.     move            y:(r2)+,y1    ;get dB-offset
  312.     cmp    a,b        #-0.6228,a    ;log of zero:  -120 dB
  313.     jle    fft_logzero
  314.     movec    m7,r7            ;#$ffff
  315.     rep    #23            ;normalize (0.5 ... 1.0)
  316.     norm    r7,b            ;shift left, decrement r7 if needed
  317.     move    b,x0
  318.     mpyr    x0,x0,a        y:(r2)+,y0    ;x**2, get c1
  319.     mpy    x0,y0,a    a,x1    y:(r2)+,y0    ;c1*x, mv x**2, get c2
  320.     mac    x1,y0,a        y:(r2)+,y0    ;c2* x**2, get c0
  321.     add    y0,a                ;add in c0
  322.     asl    a                ;multiply by 4
  323.     asl    a
  324.     asl    a                ;shift out sign bit
  325.     move    r7,a2            ;new sign = characteristic
  326.     rep    #6
  327.     asr    a            ;divide by 32, create sign bit
  328. fft_logzero
  329.     add    y1,a    #<0,y0        ;+ dB-offset
  330.     tmi    y0,a
  331.     move    #0.8027486,y0        ;1/0.31143*0.25  (display 60 dB range)
  332.     move    a,x0
  333.     mpy    x0,y0,a
  334.     asl    a
  335.     asl    a            ;*4
  336.     move    a,x:(r1)+        ;save result
  337.  
  338. fft_magnitude
  339.  
  340.     move    #<0,x0
  341.     move    x0,x:x_busy        ;clear BUSY
  342.     rts
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  
  351. ;***************************************************************************
  352.  
  353.     org    x:
  354. ;->xdata+0
  355. x_busy        ds    1
  356. ;->xdata+1
  357. x_avg_l        ds    1
  358. ;->xdata+2
  359. x_avg_r        ds    1
  360.  
  361.  
  362.  
  363.  
  364.     org    y:
  365. ;->ydata+0
  366. y_mod_count    dc    -1
  367.  
  368. y_bufpt        dc    l_buf24
  369. y_plot_type    dc    0
  370. y_plot_jmp    dc    0,plot_1,plot_2
  371.  
  372. y_fft_log_coef    dc    0.41524    ;80dB        ;0.31143==60dB
  373.         dc    0.9981958, -0.3372223, -0.6626105
  374.  
  375. ;Blackman-Harris window
  376. y_fft_window
  377.  dc 0.000000,0.000003,0.000013,0.000030,0.000054,0.000084,0.000122,0.000166
  378.  dc 0.000217,0.000274,0.000339,0.000410,0.000488,0.000573,0.000665,0.000764
  379.  dc 0.000870,0.000983,0.001102,0.001229,0.001362,0.001503,0.001650,0.001805
  380.  dc 0.001966,0.002135,0.002311,0.002494,0.002684,0.002882,0.003086,0.003298
  381.  dc 0.003517,0.003744,0.003978,0.004219,0.004468,0.004724,0.004987,0.005259
  382.  dc 0.005538,0.005824,0.006118,0.006420,0.006730,0.007047,0.007373,0.007706
  383.  dc 0.008047,0.008396,0.008753,0.009119,0.009492,0.009874,0.010264,0.010662
  384.  dc 0.011068,0.011483,0.011907,0.012339,0.012779,0.013228,0.013686,0.014153
  385.  dc 0.014628,0.015113,0.015606,0.016108,0.016619,0.017140,0.017669,0.018208
  386.  dc 0.018756,0.019314,0.019881,0.020457,0.021043,0.021639,0.022244,0.022860
  387.  dc 0.023485,0.024119,0.024764,0.025419,0.026084,0.026759,0.027445,0.028141
  388.  dc 0.028847,0.029564,0.030291,0.031029,0.031777,0.032536,0.033307,0.034087
  389.  dc 0.034879,0.035682,0.036497,0.037322,0.038158,0.039006,0.039865,0.040736
  390.  dc 0.041619,0.042512,0.043418,0.044336,0.045265,0.046206,0.047159,0.048124
  391.  dc 0.049102,0.050091,0.051093,0.052107,0.053134,0.054173,0.055224,0.056288
  392.  dc 0.057365,0.058455,0.059558,0.060673,0.061801,0.062943,0.064097,0.065265
  393.  dc 0.066446,0.067640,0.068848,0.070069,0.071304,0.072552,0.073814,0.075089
  394.  dc 0.076379,0.077682,0.078999,0.080330,0.081675,0.083034,0.084407,0.085794
  395.  dc 0.087196,0.088612,0.090042,0.091486,0.092945,0.094419,0.095907,0.097410
  396.  dc 0.098927,0.100460,0.102007,0.103568,0.105145,0.106736,0.108342,0.109964
  397.  dc 0.111600,0.113252,0.114918,0.116600,0.118297,0.120009,0.121736,0.123479
  398.  dc 0.125237,0.127010,0.128799,0.130603,0.132423,0.134258,0.136108,0.137974
  399.  dc 0.139856,0.141753,0.143666,0.145594,0.147538,0.149498,0.151473,0.153464
  400.  dc 0.155471,0.157493,0.159531,0.161585,0.163655,0.165740,0.167841,0.169957
  401.  dc 0.172090,0.174238,0.176402,0.178581,0.180776,0.182987,0.185214,0.187456
  402.  dc 0.189714,0.191988,0.194277,0.196582,0.198903,0.201239,0.203590,0.205957
  403.  dc 0.208340,0.210738,0.213152,0.215581,0.218025,0.220485,0.222960,0.225451
  404.  dc 0.227956,0.230477,0.233013,0.235564,0.238130,0.240711,0.243308,0.245919
  405.  dc 0.248544,0.251185,0.253841,0.256511,0.259196,0.261895,0.264609,0.267337
  406.  dc 0.270080,0.272836,0.275608,0.278393,0.281192,0.284005,0.286833,0.289674
  407.  dc 0.292529,0.295397,0.298279,0.301175,0.304084,0.307006,0.309941,0.312890
  408.  dc 0.315851,0.318826,0.321813,0.324813,0.327826,0.330851,0.333888,0.336938
  409.  dc 0.340000,0.343074,0.346160,0.349258,0.352367,0.355488,0.358620,0.361764
  410.  dc 0.364919,0.368085,0.371262,0.374450,0.377648,0.380857,0.384076,0.387306
  411.  dc 0.390546,0.393795,0.397055,0.400324,0.403603,0.406891,0.410188,0.413495
  412.  dc 0.416810,0.420134,0.423467,0.426808,0.430157,0.433515,0.436881,0.440254
  413.  dc 0.443635,0.447023,0.450419,0.453822,0.457231,0.460648,0.464071,0.467501
  414.  dc 0.470936,0.474378,0.477826,0.481279,0.484738,0.488202,0.491672,0.495146
  415.  dc 0.498625,0.502108,0.505596,0.509088,0.512584,0.516084,0.519587,0.523094
  416.  dc 0.526604,0.530117,0.533633,0.537151,0.540671,0.544194,0.547719,0.551245
  417.  dc 0.554773,0.558302,0.561833,0.565364,0.568896,0.572428,0.575961,0.579494
  418.  dc 0.583026,0.586558,0.590090,0.593620,0.597150,0.600678,0.604205,0.607730
  419.  dc 0.611253,0.614774,0.618292,0.621808,0.625321,0.628831,0.632337,0.635840
  420.  dc 0.639340,0.642835,0.646326,0.649813,0.653294,0.656771,0.660243,0.663710
  421.  dc 0.667170,0.670625,0.674074,0.677517,0.680953,0.684382,0.687804,0.691219
  422.  dc 0.694627,0.698027,0.701419,0.704802,0.708177,0.711544,0.714902,0.718250
  423.  dc 0.721589,0.724919,0.728239,0.731549,0.734848,0.738137,0.741415,0.744682
  424.  dc 0.747938,0.751183,0.754415,0.757636,0.760845,0.764041,0.767225,0.770396
  425.  dc 0.773553,0.776698,0.779829,0.782946,0.786049,0.789138,0.792212,0.795272
  426.  dc 0.798317,0.801347,0.804361,0.807360,0.810343,0.813310,0.816261,0.819195
  427.  dc 0.822112,0.825013,0.827897,0.830763,0.833612,0.836443,0.839255,0.842050
  428.  dc 0.844827,0.847584,0.850323,0.853043,0.855744,0.858425,0.861086,0.863728
  429.  dc 0.866349,0.868951,0.871532,0.874092,0.876631,0.879149,0.881647,0.884122
  430.  dc 0.886576,0.889008,0.891418,0.893806,0.896172,0.898515,0.900835,0.903132
  431.  dc 0.905406,0.907657,0.909884,0.912088,0.914268,0.916424,0.918556,0.920663
  432.  dc 0.922746,0.924804,0.926838,0.928846,0.930830,0.932788,0.934720,0.936627
  433.  dc 0.938508,0.940364,0.942193,0.943996,0.945772,0.947523,0.949246,0.950943
  434.  dc 0.952613,0.954256,0.955871,0.957460,0.959021,0.960554,0.962060,0.963538
  435.  dc 0.964988,0.966410,0.967804,0.969169,0.970506,0.971815,0.973095,0.974347
  436.  dc 0.975569,0.976763,0.977928,0.979064,0.980170,0.981248,0.982295,0.983314
  437.  dc 0.984303,0.985262,0.986192,0.987092,0.987962,0.988803,0.989613,0.990393
  438.  dc 0.991143,0.991864,0.992553,0.993213,0.993842,0.994441,0.995010,0.995548
  439.  dc 0.996055,0.996532,0.996979,0.997394,0.997779,0.998134,0.998457,0.998750
  440.  dc 0.999013,0.999244,0.999444,0.999614,0.999753,0.999861,0.999938,0.999985
  441.  dc 0.999999,0.999985,0.999938,0.999861,0.999753,0.999614,0.999444,0.999244
  442.  dc 0.999013,0.998750,0.998457,0.998134,0.997779,0.997394,0.996979,0.996532
  443.  dc 0.996055,0.995548,0.995010,0.994441,0.993842,0.993213,0.992553,0.991864
  444.  dc 0.991143,0.990393,0.989613,0.988803,0.987962,0.987092,0.986192,0.985262
  445.  dc 0.984303,0.983314,0.982295,0.981248,0.980170,0.979064,0.977928,0.976763
  446.  dc 0.975569,0.974347,0.973095,0.971815,0.970506,0.969169,0.967804,0.966410
  447.  dc 0.964988,0.963538,0.962060,0.960554,0.959021,0.957460,0.955871,0.954256
  448.  dc 0.952613,0.950943,0.949246,0.947523,0.945772,0.943996,0.942193,0.940364
  449.  dc 0.938508,0.936627,0.934720,0.932788,0.930830,0.928846,0.926838,0.924804
  450.  dc 0.922746,0.920663,0.918556,0.916424,0.914268,0.912088,0.909884,0.907657
  451.  dc 0.905406,0.903132,0.900835,0.898515,0.896172,0.893806,0.891418,0.889008
  452.  dc 0.886576,0.884122,0.881647,0.879149,0.876631,0.874092,0.871532,0.868951
  453.  dc 0.866349,0.863728,0.861086,0.858425,0.855744,0.853043,0.850323,0.847584
  454.  dc 0.844827,0.842050,0.839255,0.836443,0.833612,0.830763,0.827897,0.825013
  455.  dc 0.822112,0.819195,0.816261,0.813310,0.810343,0.807360,0.804361,0.801347
  456.  dc 0.798317,0.795272,0.792212,0.789138,0.786049,0.782946,0.779829,0.776698
  457.  dc 0.773553,0.770396,0.767225,0.764041,0.760845,0.757636,0.754415,0.751183
  458.  dc 0.747938,0.744682,0.741415,0.738137,0.734848,0.731549,0.728239,0.724919
  459.  dc 0.721589,0.718250,0.714902,0.711544,0.708177,0.704802,0.701419,0.698027
  460.  dc 0.694627,0.691219,0.687804,0.684382,0.680953,0.677517,0.674074,0.670625
  461.  dc 0.667170,0.663710,0.660243,0.656771,0.653294,0.649813,0.646326,0.642835
  462.  dc 0.639340,0.635840,0.632337,0.628831,0.625321,0.621808,0.618292,0.614774
  463.  dc 0.611253,0.607730,0.604205,0.600678,0.597150,0.593620,0.590090,0.586558
  464.  dc 0.583026,0.579494,0.575961,0.572428,0.568896,0.565364,0.561833,0.558302
  465.  dc 0.554773,0.551245,0.547719,0.544194,0.540671,0.537151,0.533633,0.530117
  466.  dc 0.526604,0.523094,0.519587,0.516084,0.512584,0.509088,0.505596,0.502108
  467.  dc 0.498625,0.495146,0.491672,0.488202,0.484738,0.481279,0.477826,0.474378
  468.  dc 0.470936,0.467501,0.464071,0.460648,0.457231,0.453822,0.450419,0.447023
  469.  dc 0.443635,0.440254,0.436881,0.433515,0.430157,0.426808,0.423467,0.420134
  470.  dc 0.416810,0.413495,0.410188,0.406891,0.403603,0.400324,0.397055,0.393795
  471.  dc 0.390546,0.387306,0.384076,0.380857,0.377648,0.374450,0.371262,0.368085
  472.  dc 0.364919,0.361764,0.358620,0.355488,0.352367,0.349258,0.346160,0.343074
  473.  dc 0.340000,0.336938,0.333888,0.330851,0.327826,0.324813,0.321813,0.318826
  474.  dc 0.315851,0.312890,0.309941,0.307006,0.304084,0.301175,0.298279,0.295397
  475.  dc 0.292529,0.289674,0.286833,0.284005,0.281192,0.278393,0.275608,0.272836
  476.  dc 0.270080,0.267337,0.264609,0.261895,0.259196,0.256511,0.253841,0.251185
  477.  dc 0.248544,0.245919,0.243308,0.240711,0.238130,0.235564,0.233013,0.230477
  478.  dc 0.227956,0.225451,0.222960,0.220485,0.218025,0.215581,0.213152,0.210738
  479.  dc 0.208340,0.205957,0.203590,0.201239,0.198903,0.196582,0.194277,0.191988
  480.  dc 0.189714,0.187456,0.185214,0.182987,0.180776,0.178581,0.176402,0.174238
  481.  dc 0.172090,0.169957,0.167841,0.165740,0.163655,0.161585,0.159531,0.157493
  482.  dc 0.155471,0.153464,0.151473,0.149498,0.147538,0.145594,0.143666,0.141753
  483.  dc 0.139856,0.137974,0.136108,0.134258,0.132423,0.130603,0.128799,0.127010
  484.  dc 0.125237,0.123479,0.121736,0.120009,0.118297,0.116600,0.114918,0.113252
  485.  dc 0.111600,0.109964,0.108342,0.106736,0.105145,0.103568,0.102007,0.100460
  486.  dc 0.098927,0.097410,0.095907,0.094419,0.092945,0.091486,0.090042,0.088612
  487.  dc 0.087196,0.085794,0.084407,0.083034,0.081675,0.080330,0.078999,0.077682
  488.  dc 0.076379,0.075089,0.073814,0.072552,0.071304,0.070069,0.068848,0.067640
  489.  dc 0.066446,0.065265,0.064097,0.062943,0.061801,0.060673,0.059558,0.058455
  490.  dc 0.057365,0.056288,0.055224,0.054173,0.053134,0.052107,0.051093,0.050091
  491.  dc 0.049102,0.048124,0.047159,0.046206,0.045265,0.044336,0.043418,0.042512
  492.  dc 0.041619,0.040736,0.039865,0.039006,0.038158,0.037322,0.036497,0.035682
  493.  dc 0.034879,0.034087,0.033307,0.032536,0.031777,0.031029,0.030291,0.029564
  494.  dc 0.028847,0.028141,0.027445,0.026759,0.026084,0.025419,0.024764,0.024119
  495.  dc 0.023485,0.022860,0.022244,0.021639,0.021043,0.020457,0.019881,0.019314
  496.  dc 0.018756,0.018208,0.017669,0.017140,0.016619,0.016108,0.015606,0.015113
  497.  dc 0.014628,0.014153,0.013686,0.013228,0.012779,0.012339,0.011907,0.011483
  498.  dc 0.011068,0.010662,0.010264,0.009874,0.009492,0.009119,0.008753,0.008396
  499.  dc 0.008047,0.007706,0.007373,0.007047,0.006730,0.006420,0.006118,0.005824
  500.  dc 0.005538,0.005259,0.004987,0.004724,0.004468,0.004219,0.003978,0.003744
  501.  dc 0.003517,0.003298,0.003086,0.002882,0.002684,0.002494,0.002311,0.002135
  502.  dc 0.001966,0.001805,0.001650,0.001503,0.001362,0.001229,0.001102,0.000983
  503.  dc 0.000870,0.000764,0.000665,0.000573,0.000488,0.000410,0.000339,0.000274
  504.  dc 0.000217,0.000166,0.000122,0.000084,0.000054,0.000030,0.000013,0.000003
  505.  
  506.  
  507.  
  508.  
  509.     org    l:
  510.     align    2048
  511. l_buf24        ds    2048        ;ALIGN 2048 required!
  512.  
  513. ;->ldata+2048
  514. l_fft_out
  515. l_buf8        ds    1024        ;ALIGN 1024 required!
  516.  
  517. l_fft_in    ds    1024        ;ALIGN 1024 required!
  518.  
  519.  
  520. l_fft_sincos    ;x:-cos        ;y:-sin    ;ALIGN 512  required!
  521.     dcl    -1.000000,    0.000000    ;0
  522.     dcl    -0.999981,    -0.006135    ;1
  523.     dcl    -0.999925,    -0.012271    ;2
  524.     dcl    -0.999831,    -0.018406    ;3
  525.     dcl    -0.999699,    -0.024541    ;4
  526.     dcl    -0.999529,    -0.030674    ;5
  527.     dcl    -0.999322,    -0.036807    ;6
  528.     dcl    -0.999078,    -0.042938    ;7
  529.     dcl    -0.998795,    -0.049067    ;8
  530.     dcl    -0.998476,    -0.055195    ;9
  531.     dcl    -0.998118,    -0.061320    ;10
  532.     dcl    -0.997723,    -0.067443    ;11
  533.     dcl    -0.997290,    -0.073564    ;12
  534.     dcl    -0.996820,    -0.079682    ;13
  535.     dcl    -0.996313,    -0.085797    ;14
  536.     dcl    -0.995767,    -0.091909    ;15
  537.     dcl    -0.995185,    -0.098017    ;16
  538.     dcl    -0.994565,    -0.104122    ;17
  539.     dcl    -0.993907,    -0.110222    ;18
  540.     dcl    -0.993212,    -0.116319    ;19
  541.     dcl    -0.992480,    -0.122411    ;20
  542.     dcl    -0.991710,    -0.128498    ;21
  543.     dcl    -0.990903,    -0.134581    ;22
  544.     dcl    -0.990058,    -0.140658    ;23
  545.     dcl    -0.989177,    -0.146730    ;24
  546.     dcl    -0.988258,    -0.152797    ;25
  547.     dcl    -0.987301,    -0.158858    ;26
  548.     dcl    -0.986308,    -0.164913    ;27
  549.     dcl    -0.985278,    -0.170962    ;28
  550.     dcl    -0.984210,    -0.177004    ;29
  551.     dcl    -0.983105,    -0.183040    ;30
  552.     dcl    -0.981964,    -0.189069    ;31
  553.     dcl    -0.980785,    -0.195090    ;32
  554.     dcl    -0.979570,    -0.201105    ;33
  555.     dcl    -0.978317,    -0.207111    ;34
  556.     dcl    -0.977028,    -0.213110    ;35
  557.     dcl    -0.975702,    -0.219101    ;36
  558.     dcl    -0.974339,    -0.225084    ;37
  559.     dcl    -0.972940,    -0.231058    ;38
  560.     dcl    -0.971504,    -0.237024    ;39
  561.     dcl    -0.970031,    -0.242980    ;40
  562.     dcl    -0.968522,    -0.248928    ;41
  563.     dcl    -0.966976,    -0.254866    ;42
  564.     dcl    -0.965394,    -0.260794    ;43
  565.     dcl    -0.963776,    -0.266713    ;44
  566.     dcl    -0.962121,    -0.272621    ;45
  567.     dcl    -0.960431,    -0.278520    ;46
  568.     dcl    -0.958703,    -0.284408    ;47
  569.     dcl    -0.956940,    -0.290285    ;48
  570.     dcl    -0.955141,    -0.296151    ;49
  571.     dcl    -0.953306,    -0.302006    ;50
  572.     dcl    -0.951435,    -0.307850    ;51
  573.     dcl    -0.949528,    -0.313682    ;52
  574.     dcl    -0.947586,    -0.319502    ;53
  575.     dcl    -0.945607,    -0.325310    ;54
  576.     dcl    -0.943593,    -0.331106    ;55
  577.     dcl    -0.941544,    -0.336890    ;56
  578.     dcl    -0.939459,    -0.342661    ;57
  579.     dcl    -0.937339,    -0.348419    ;58
  580.     dcl    -0.935184,    -0.354164    ;59
  581.     dcl    -0.932993,    -0.359895    ;60
  582.     dcl    -0.930767,    -0.365613    ;61
  583.     dcl    -0.928506,    -0.371317    ;62
  584.     dcl    -0.926210,    -0.377007    ;63
  585.     dcl    -0.923880,    -0.382683    ;64
  586.     dcl    -0.921514,    -0.388345    ;65
  587.     dcl    -0.919114,    -0.393992    ;66
  588.     dcl    -0.916679,    -0.399624    ;67
  589.     dcl    -0.914210,    -0.405241    ;68
  590.     dcl    -0.911706,    -0.410843    ;69
  591.     dcl    -0.909168,    -0.416430    ;70
  592.     dcl    -0.906596,    -0.422000    ;71
  593.     dcl    -0.903989,    -0.427555    ;72
  594.     dcl    -0.901349,    -0.433094    ;73
  595.     dcl    -0.898674,    -0.438616    ;74
  596.     dcl    -0.895966,    -0.444122    ;75
  597.     dcl    -0.893224,    -0.449611    ;76
  598.     dcl    -0.890449,    -0.455084    ;77
  599.     dcl    -0.887640,    -0.460539    ;78
  600.     dcl    -0.884797,    -0.465976    ;79
  601.     dcl    -0.881921,    -0.471397    ;80
  602.     dcl    -0.879012,    -0.476799    ;81
  603.     dcl    -0.876070,    -0.482184    ;82
  604.     dcl    -0.873095,    -0.487550    ;83
  605.     dcl    -0.870087,    -0.492898    ;84
  606.     dcl    -0.867046,    -0.498228    ;85
  607.     dcl    -0.863973,    -0.503538    ;86
  608.     dcl    -0.860867,    -0.508830    ;87
  609.     dcl    -0.857729,    -0.514103    ;88
  610.     dcl    -0.854558,    -0.519356    ;89
  611.     dcl    -0.851355,    -0.524590    ;90
  612.     dcl    -0.848120,    -0.529804    ;91
  613.     dcl    -0.844854,    -0.534998    ;92
  614.     dcl    -0.841555,    -0.540171    ;93
  615.     dcl    -0.838225,    -0.545325    ;94
  616.     dcl    -0.834863,    -0.550458    ;95
  617.     dcl    -0.831470,    -0.555570    ;96
  618.     dcl    -0.828045,    -0.560662    ;97
  619.     dcl    -0.824589,    -0.565732    ;98
  620.     dcl    -0.821103,    -0.570781    ;99
  621.     dcl    -0.817585,    -0.575808    ;100
  622.     dcl    -0.814036,    -0.580814    ;101
  623.     dcl    -0.810457,    -0.585798    ;102
  624.     dcl    -0.806848,    -0.590760    ;103
  625.     dcl    -0.803208,    -0.595699    ;104
  626.     dcl    -0.799537,    -0.600616    ;105
  627.     dcl    -0.795837,    -0.605511    ;106
  628.     dcl    -0.792107,    -0.610383    ;107
  629.     dcl    -0.788346,    -0.615232    ;108
  630.     dcl    -0.784557,    -0.620057    ;109
  631.     dcl    -0.780737,    -0.624859    ;110
  632.     dcl    -0.776888,    -0.629638    ;111
  633.     dcl    -0.773010,    -0.634393    ;112
  634.     dcl    -0.769103,    -0.639124    ;113
  635.     dcl    -0.765167,    -0.643832    ;114
  636.     dcl    -0.761202,    -0.648514    ;115
  637.     dcl    -0.757209,    -0.653173    ;116
  638.     dcl    -0.753187,    -0.657807    ;117
  639.     dcl    -0.749136,    -0.662416    ;118
  640.     dcl    -0.745058,    -0.667000    ;119
  641.     dcl    -0.740951,    -0.671559    ;120
  642.     dcl    -0.736817,    -0.676093    ;121
  643.     dcl    -0.732654,    -0.680601    ;122
  644.     dcl    -0.728464,    -0.685084    ;123
  645.     dcl    -0.724247,    -0.689541    ;124
  646.     dcl    -0.720003,    -0.693971    ;125
  647.     dcl    -0.715731,    -0.698376    ;126
  648.     dcl    -0.711432,    -0.702755    ;127
  649.     dcl    -0.707107,    -0.707107    ;128
  650.     dcl    -0.702755,    -0.711432    ;129
  651.     dcl    -0.698376,    -0.715731    ;130
  652.     dcl    -0.693971,    -0.720003    ;131
  653.     dcl    -0.689541,    -0.724247    ;132
  654.     dcl    -0.685084,    -0.728464    ;133
  655.     dcl    -0.680601,    -0.732654    ;134
  656.     dcl    -0.676093,    -0.736817    ;135
  657.     dcl    -0.671559,    -0.740951    ;136
  658.     dcl    -0.667000,    -0.745058    ;137
  659.     dcl    -0.662416,    -0.749136    ;138
  660.     dcl    -0.657807,    -0.753187    ;139
  661.     dcl    -0.653173,    -0.757209    ;140
  662.     dcl    -0.648514,    -0.761202    ;141
  663.     dcl    -0.643832,    -0.765167    ;142
  664.     dcl    -0.639124,    -0.769103    ;143
  665.     dcl    -0.634393,    -0.773010    ;144
  666.     dcl    -0.629638,    -0.776888    ;145
  667.     dcl    -0.624859,    -0.780737    ;146
  668.     dcl    -0.620057,    -0.784557    ;147
  669.     dcl    -0.615232,    -0.788346    ;148
  670.     dcl    -0.610383,    -0.792107    ;149
  671.     dcl    -0.605511,    -0.795837    ;150
  672.     dcl    -0.600616,    -0.799537    ;151
  673.     dcl    -0.595699,    -0.803208    ;152
  674.     dcl    -0.590760,    -0.806848    ;153
  675.     dcl    -0.585798,    -0.810457    ;154
  676.     dcl    -0.580814,    -0.814036    ;155
  677.     dcl    -0.575808,    -0.817585    ;156
  678.     dcl    -0.570781,    -0.821103    ;157
  679.     dcl    -0.565732,    -0.824589    ;158
  680.     dcl    -0.560662,    -0.828045    ;159
  681.     dcl    -0.555570,    -0.831470    ;160
  682.     dcl    -0.550458,    -0.834863    ;161
  683.     dcl    -0.545325,    -0.838225    ;162
  684.     dcl    -0.540171,    -0.841555    ;163
  685.     dcl    -0.534998,    -0.844854    ;164
  686.     dcl    -0.529804,    -0.848120    ;165
  687.     dcl    -0.524590,    -0.851355    ;166
  688.     dcl    -0.519356,    -0.854558    ;167
  689.     dcl    -0.514103,    -0.857729    ;168
  690.     dcl    -0.508830,    -0.860867    ;169
  691.     dcl    -0.503538,    -0.863973    ;170
  692.     dcl    -0.498228,    -0.867046    ;171
  693.     dcl    -0.492898,    -0.870087    ;172
  694.     dcl    -0.487550,    -0.873095    ;173
  695.     dcl    -0.482184,    -0.876070    ;174
  696.     dcl    -0.476799,    -0.879012    ;175
  697.     dcl    -0.471397,    -0.881921    ;176
  698.     dcl    -0.465976,    -0.884797    ;177
  699.     dcl    -0.460539,    -0.887640    ;178
  700.     dcl    -0.455084,    -0.890449    ;179
  701.     dcl    -0.449611,    -0.893224    ;180
  702.     dcl    -0.444122,    -0.895966    ;181
  703.     dcl    -0.438616,    -0.898674    ;182
  704.     dcl    -0.433094,    -0.901349    ;183
  705.     dcl    -0.427555,    -0.903989    ;184
  706.     dcl    -0.422000,    -0.906596    ;185
  707.     dcl    -0.416430,    -0.909168    ;186
  708.     dcl    -0.410843,    -0.911706    ;187
  709.     dcl    -0.405241,    -0.914210    ;188
  710.     dcl    -0.399624,    -0.916679    ;189
  711.     dcl    -0.393992,    -0.919114    ;190
  712.     dcl    -0.388345,    -0.921514    ;191
  713.     dcl    -0.382683,    -0.923880    ;192
  714.     dcl    -0.377007,    -0.926210    ;193
  715.     dcl    -0.371317,    -0.928506    ;194
  716.     dcl    -0.365613,    -0.930767    ;195
  717.     dcl    -0.359895,    -0.932993    ;196
  718.     dcl    -0.354164,    -0.935184    ;197
  719.     dcl    -0.348419,    -0.937339    ;198
  720.     dcl    -0.342661,    -0.939459    ;199
  721.     dcl    -0.336890,    -0.941544    ;200
  722.     dcl    -0.331106,    -0.943593    ;201
  723.     dcl    -0.325310,    -0.945607    ;202
  724.     dcl    -0.319502,    -0.947586    ;203
  725.     dcl    -0.313682,    -0.949528    ;204
  726.     dcl    -0.307850,    -0.951435    ;205
  727.     dcl    -0.302006,    -0.953306    ;206
  728.     dcl    -0.296151,    -0.955141    ;207
  729.     dcl    -0.290285,    -0.956940    ;208
  730.     dcl    -0.284408,    -0.958703    ;209
  731.     dcl    -0.278520,    -0.960431    ;210
  732.     dcl    -0.272621,    -0.962121    ;211
  733.     dcl    -0.266713,    -0.963776    ;212
  734.     dcl    -0.260794,    -0.965394    ;213
  735.     dcl    -0.254866,    -0.966976    ;214
  736.     dcl    -0.248928,    -0.968522    ;215
  737.     dcl    -0.242980,    -0.970031    ;216
  738.     dcl    -0.237024,    -0.971504    ;217
  739.     dcl    -0.231058,    -0.972940    ;218
  740.     dcl    -0.225084,    -0.974339    ;219
  741.     dcl    -0.219101,    -0.975702    ;220
  742.     dcl    -0.213110,    -0.977028    ;221
  743.     dcl    -0.207111,    -0.978317    ;222
  744.     dcl    -0.201105,    -0.979570    ;223
  745.     dcl    -0.195090,    -0.980785    ;224
  746.     dcl    -0.189069,    -0.981964    ;225
  747.     dcl    -0.183040,    -0.983105    ;226
  748.     dcl    -0.177004,    -0.984210    ;227
  749.     dcl    -0.170962,    -0.985278    ;228
  750.     dcl    -0.164913,    -0.986308    ;229
  751.     dcl    -0.158858,    -0.987301    ;230
  752.     dcl    -0.152797,    -0.988258    ;231
  753.     dcl    -0.146730,    -0.989177    ;232
  754.     dcl    -0.140658,    -0.990058    ;233
  755.     dcl    -0.134581,    -0.990903    ;234
  756.     dcl    -0.128498,    -0.991710    ;235
  757.     dcl    -0.122411,    -0.992480    ;236
  758.     dcl    -0.116319,    -0.993212    ;237
  759.     dcl    -0.110222,    -0.993907    ;238
  760.     dcl    -0.104122,    -0.994565    ;239
  761.     dcl    -0.098017,    -0.995185    ;240
  762.     dcl    -0.091909,    -0.995767    ;241
  763.     dcl    -0.085797,    -0.996313    ;242
  764.     dcl    -0.079682,    -0.996820    ;243
  765.     dcl    -0.073564,    -0.997290    ;244
  766.     dcl    -0.067443,    -0.997723    ;245
  767.     dcl    -0.061320,    -0.998118    ;246
  768.     dcl    -0.055195,    -0.998476    ;247
  769.     dcl    -0.049067,    -0.998795    ;248
  770.     dcl    -0.042938,    -0.999078    ;249
  771.     dcl    -0.036807,    -0.999322    ;250
  772.     dcl    -0.030674,    -0.999529    ;251
  773.     dcl    -0.024541,    -0.999699    ;252
  774.     dcl    -0.018406,    -0.999831    ;253
  775.     dcl    -0.012271,    -0.999925    ;254
  776.     dcl    -0.006135,    -0.999981    ;255
  777.     dcl    0.000000,    -1.000000    ;256
  778.     dcl    0.006135,    -0.999981    ;257
  779.     dcl    0.012271,    -0.999925    ;258
  780.     dcl    0.018406,    -0.999831    ;259
  781.     dcl    0.024541,    -0.999699    ;260
  782.     dcl    0.030674,    -0.999529    ;261
  783.     dcl    0.036807,    -0.999322    ;262
  784.     dcl    0.042938,    -0.999078    ;263
  785.     dcl    0.049067,    -0.998795    ;264
  786.     dcl    0.055195,    -0.998476    ;265
  787.     dcl    0.061320,    -0.998118    ;266
  788.     dcl    0.067443,    -0.997723    ;267
  789.     dcl    0.073564,    -0.997290    ;268
  790.     dcl    0.079682,    -0.996820    ;269
  791.     dcl    0.085797,    -0.996313    ;270
  792.     dcl    0.091909,    -0.995767    ;271
  793.     dcl    0.098017,    -0.995185    ;272
  794.     dcl    0.104122,    -0.994565    ;273
  795.     dcl    0.110222,    -0.993907    ;274
  796.     dcl    0.116319,    -0.993212    ;275
  797.     dcl    0.122411,    -0.992480    ;276
  798.     dcl    0.128498,    -0.991710    ;277
  799.     dcl    0.134581,    -0.990903    ;278
  800.     dcl    0.140658,    -0.990058    ;279
  801.     dcl    0.146730,    -0.989177    ;280
  802.     dcl    0.152797,    -0.988258    ;281
  803.     dcl    0.158858,    -0.987301    ;282
  804.     dcl    0.164913,    -0.986308    ;283
  805.     dcl    0.170962,    -0.985278    ;284
  806.     dcl    0.177004,    -0.984210    ;285
  807.     dcl    0.183040,    -0.983105    ;286
  808.     dcl    0.189069,    -0.981964    ;287
  809.     dcl    0.195090,    -0.980785    ;288
  810.     dcl    0.201105,    -0.979570    ;289
  811.     dcl    0.207111,    -0.978317    ;290
  812.     dcl    0.213110,    -0.977028    ;291
  813.     dcl    0.219101,    -0.975702    ;292
  814.     dcl    0.225084,    -0.974339    ;293
  815.     dcl    0.231058,    -0.972940    ;294
  816.     dcl    0.237024,    -0.971504    ;295
  817.     dcl    0.242980,    -0.970031    ;296
  818.     dcl    0.248928,    -0.968522    ;297
  819.     dcl    0.254866,    -0.966976    ;298
  820.     dcl    0.260794,    -0.965394    ;299
  821.     dcl    0.266713,    -0.963776    ;300
  822.     dcl    0.272621,    -0.962121    ;301
  823.     dcl    0.278520,    -0.960431    ;302
  824.     dcl    0.284408,    -0.958703    ;303
  825.     dcl    0.290285,    -0.956940    ;304
  826.     dcl    0.296151,    -0.955141    ;305
  827.     dcl    0.302006,    -0.953306    ;306
  828.     dcl    0.307850,    -0.951435    ;307
  829.     dcl    0.313682,    -0.949528    ;308
  830.     dcl    0.319502,    -0.947586    ;309
  831.     dcl    0.325310,    -0.945607    ;310
  832.     dcl    0.331106,    -0.943593    ;311
  833.     dcl    0.336890,    -0.941544    ;312
  834.     dcl    0.342661,    -0.939459    ;313
  835.     dcl    0.348419,    -0.937339    ;314
  836.     dcl    0.354164,    -0.935184    ;315
  837.     dcl    0.359895,    -0.932993    ;316
  838.     dcl    0.365613,    -0.930767    ;317
  839.     dcl    0.371317,    -0.928506    ;318
  840.     dcl    0.377007,    -0.926210    ;319
  841.     dcl    0.382683,    -0.923880    ;320
  842.     dcl    0.388345,    -0.921514    ;321
  843.     dcl    0.393992,    -0.919114    ;322
  844.     dcl    0.399624,    -0.916679    ;323
  845.     dcl    0.405241,    -0.914210    ;324
  846.     dcl    0.410843,    -0.911706    ;325
  847.     dcl    0.416430,    -0.909168    ;326
  848.     dcl    0.422000,    -0.906596    ;327
  849.     dcl    0.427555,    -0.903989    ;328
  850.     dcl    0.433094,    -0.901349    ;329
  851.     dcl    0.438616,    -0.898674    ;330
  852.     dcl    0.444122,    -0.895966    ;331
  853.     dcl    0.449611,    -0.893224    ;332
  854.     dcl    0.455084,    -0.890449    ;333
  855.     dcl    0.460539,    -0.887640    ;334
  856.     dcl    0.465976,    -0.884797    ;335
  857.     dcl    0.471397,    -0.881921    ;336
  858.     dcl    0.476799,    -0.879012    ;337
  859.     dcl    0.482184,    -0.876070    ;338
  860.     dcl    0.487550,    -0.873095    ;339
  861.     dcl    0.492898,    -0.870087    ;340
  862.     dcl    0.498228,    -0.867046    ;341
  863.     dcl    0.503538,    -0.863973    ;342
  864.     dcl    0.508830,    -0.860867    ;343
  865.     dcl    0.514103,    -0.857729    ;344
  866.     dcl    0.519356,    -0.854558    ;345
  867.     dcl    0.524590,    -0.851355    ;346
  868.     dcl    0.529804,    -0.848120    ;347
  869.     dcl    0.534998,    -0.844854    ;348
  870.     dcl    0.540171,    -0.841555    ;349
  871.     dcl    0.545325,    -0.838225    ;350
  872.     dcl    0.550458,    -0.834863    ;351
  873.     dcl    0.555570,    -0.831470    ;352
  874.     dcl    0.560662,    -0.828045    ;353
  875.     dcl    0.565732,    -0.824589    ;354
  876.     dcl    0.570781,    -0.821103    ;355
  877.     dcl    0.575808,    -0.817585    ;356
  878.     dcl    0.580814,    -0.814036    ;357
  879.     dcl    0.585798,    -0.810457    ;358
  880.     dcl    0.590760,    -0.806848    ;359
  881.     dcl    0.595699,    -0.803208    ;360
  882.     dcl    0.600616,    -0.799537    ;361
  883.     dcl    0.605511,    -0.795837    ;362
  884.     dcl    0.610383,    -0.792107    ;363
  885.     dcl    0.615232,    -0.788346    ;364
  886.     dcl    0.620057,    -0.784557    ;365
  887.     dcl    0.624859,    -0.780737    ;366
  888.     dcl    0.629638,    -0.776888    ;367
  889.     dcl    0.634393,    -0.773010    ;368
  890.     dcl    0.639124,    -0.769103    ;369
  891.     dcl    0.643832,    -0.765167    ;370
  892.     dcl    0.648514,    -0.761202    ;371
  893.     dcl    0.653173,    -0.757209    ;372
  894.     dcl    0.657807,    -0.753187    ;373
  895.     dcl    0.662416,    -0.749136    ;374
  896.     dcl    0.667000,    -0.745058    ;375
  897.     dcl    0.671559,    -0.740951    ;376
  898.     dcl    0.676093,    -0.736817    ;377
  899.     dcl    0.680601,    -0.732654    ;378
  900.     dcl    0.685084,    -0.728464    ;379
  901.     dcl    0.689541,    -0.724247    ;380
  902.     dcl    0.693971,    -0.720003    ;381
  903.     dcl    0.698376,    -0.715731    ;382
  904.     dcl    0.702755,    -0.711432    ;383
  905.     dcl    0.707107,    -0.707107    ;384
  906.     dcl    0.711432,    -0.702755    ;385
  907.     dcl    0.715731,    -0.698376    ;386
  908.     dcl    0.720003,    -0.693971    ;387
  909.     dcl    0.724247,    -0.689541    ;388
  910.     dcl    0.728464,    -0.685084    ;389
  911.     dcl    0.732654,    -0.680601    ;390
  912.     dcl    0.736817,    -0.676093    ;391
  913.     dcl    0.740951,    -0.671559    ;392
  914.     dcl    0.745058,    -0.667000    ;393
  915.     dcl    0.749136,    -0.662416    ;394
  916.     dcl    0.753187,    -0.657807    ;395
  917.     dcl    0.757209,    -0.653173    ;396
  918.     dcl    0.761202,    -0.648514    ;397
  919.     dcl    0.765167,    -0.643832    ;398
  920.     dcl    0.769103,    -0.639124    ;399
  921.     dcl    0.773010,    -0.634393    ;400
  922.     dcl    0.776888,    -0.629638    ;401
  923.     dcl    0.780737,    -0.624859    ;402
  924.     dcl    0.784557,    -0.620057    ;403
  925.     dcl    0.788346,    -0.615232    ;404
  926.     dcl    0.792107,    -0.610383    ;405
  927.     dcl    0.795837,    -0.605511    ;406
  928.     dcl    0.799537,    -0.600616    ;407
  929.     dcl    0.803208,    -0.595699    ;408
  930.     dcl    0.806848,    -0.590760    ;409
  931.     dcl    0.810457,    -0.585798    ;410
  932.     dcl    0.814036,    -0.580814    ;411
  933.     dcl    0.817585,    -0.575808    ;412
  934.     dcl    0.821103,    -0.570781    ;413
  935.     dcl    0.824589,    -0.565732    ;414
  936.     dcl    0.828045,    -0.560662    ;415
  937.     dcl    0.831470,    -0.555570    ;416
  938.     dcl    0.834863,    -0.550458    ;417
  939.     dcl    0.838225,    -0.545325    ;418
  940.     dcl    0.841555,    -0.540171    ;419
  941.     dcl    0.844854,    -0.534998    ;420
  942.     dcl    0.848120,    -0.529804    ;421
  943.     dcl    0.851355,    -0.524590    ;422
  944.     dcl    0.854558,    -0.519356    ;423
  945.     dcl    0.857729,    -0.514103    ;424
  946.     dcl    0.860867,    -0.508830    ;425
  947.     dcl    0.863973,    -0.503538    ;426
  948.     dcl    0.867046,    -0.498228    ;427
  949.     dcl    0.870087,    -0.492898    ;428
  950.     dcl    0.873095,    -0.487550    ;429
  951.     dcl    0.876070,    -0.482184    ;430
  952.     dcl    0.879012,    -0.476799    ;431
  953.     dcl    0.881921,    -0.471397    ;432
  954.     dcl    0.884797,    -0.465976    ;433
  955.     dcl    0.887640,    -0.460539    ;434
  956.     dcl    0.890449,    -0.455084    ;435
  957.     dcl    0.893224,    -0.449611    ;436
  958.     dcl    0.895966,    -0.444122    ;437
  959.     dcl    0.898674,    -0.438616    ;438
  960.     dcl    0.901349,    -0.433094    ;439
  961.     dcl    0.903989,    -0.427555    ;440
  962.     dcl    0.906596,    -0.422000    ;441
  963.     dcl    0.909168,    -0.416430    ;442
  964.     dcl    0.911706,    -0.410843    ;443
  965.     dcl    0.914210,    -0.405241    ;444
  966.     dcl    0.916679,    -0.399624    ;445
  967.     dcl    0.919114,    -0.393992    ;446
  968.     dcl    0.921514,    -0.388345    ;447
  969.     dcl    0.923880,    -0.382683    ;448
  970.     dcl    0.926210,    -0.377007    ;449
  971.     dcl    0.928506,    -0.371317    ;450
  972.     dcl    0.930767,    -0.365613    ;451
  973.     dcl    0.932993,    -0.359895    ;452
  974.     dcl    0.935184,    -0.354164    ;453
  975.     dcl    0.937339,    -0.348419    ;454
  976.     dcl    0.939459,    -0.342661    ;455
  977.     dcl    0.941544,    -0.336890    ;456
  978.     dcl    0.943593,    -0.331106    ;457
  979.     dcl    0.945607,    -0.325310    ;458
  980.     dcl    0.947586,    -0.319502    ;459
  981.     dcl    0.949528,    -0.313682    ;460
  982.     dcl    0.951435,    -0.307850    ;461
  983.     dcl    0.953306,    -0.302006    ;462
  984.     dcl    0.955141,    -0.296151    ;463
  985.     dcl    0.956940,    -0.290285    ;464
  986.     dcl    0.958703,    -0.284408    ;465
  987.     dcl    0.960431,    -0.278520    ;466
  988.     dcl    0.962121,    -0.272621    ;467
  989.     dcl    0.963776,    -0.266713    ;468
  990.     dcl    0.965394,    -0.260794    ;469
  991.     dcl    0.966976,    -0.254866    ;470
  992.     dcl    0.968522,    -0.248928    ;471
  993.     dcl    0.970031,    -0.242980    ;472
  994.     dcl    0.971504,    -0.237024    ;473
  995.     dcl    0.972940,    -0.231058    ;474
  996.     dcl    0.974339,    -0.225084    ;475
  997.     dcl    0.975702,    -0.219101    ;476
  998.     dcl    0.977028,    -0.213110    ;477
  999.     dcl    0.978317,    -0.207111    ;478
  1000.     dcl    0.979570,    -0.201105    ;479
  1001.     dcl    0.980785,    -0.195090    ;480
  1002.     dcl    0.981964,    -0.189069    ;481
  1003.     dcl    0.983105,    -0.183040    ;482
  1004.     dcl    0.984210,    -0.177004    ;483
  1005.     dcl    0.985278,    -0.170962    ;484
  1006.     dcl    0.986308,    -0.164913    ;485
  1007.     dcl    0.987301,    -0.158858    ;486
  1008.     dcl    0.988258,    -0.152797    ;487
  1009.     dcl    0.989177,    -0.146730    ;488
  1010.     dcl    0.990058,    -0.140658    ;489
  1011.     dcl    0.990903,    -0.134581    ;490
  1012.     dcl    0.991710,    -0.128498    ;491
  1013.     dcl    0.992480,    -0.122411    ;492
  1014.     dcl    0.993212,    -0.116319    ;493
  1015.     dcl    0.993907,    -0.110222    ;494
  1016.     dcl    0.994565,    -0.104122    ;495
  1017.     dcl    0.995185,    -0.098017    ;496
  1018.     dcl    0.995767,    -0.091909    ;497
  1019.     dcl    0.996313,    -0.085797    ;498
  1020.     dcl    0.996820,    -0.079682    ;499
  1021.     dcl    0.997290,    -0.073564    ;500
  1022.     dcl    0.997723,    -0.067443    ;501
  1023.     dcl    0.998118,    -0.061320    ;502
  1024.     dcl    0.998476,    -0.055195    ;503
  1025.     dcl    0.998795,    -0.049067    ;504
  1026.     dcl    0.999078,    -0.042938    ;505
  1027.     dcl    0.999322,    -0.036807    ;506
  1028.     dcl    0.999529,    -0.030674    ;507
  1029.     dcl    0.999699,    -0.024541    ;508
  1030.     dcl    0.999831,    -0.018406    ;509
  1031.     dcl    0.999925,    -0.012271    ;510
  1032.     dcl    0.999981,    -0.006135    ;511
  1033.  
  1034. ;the END
  1035.